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Methodology is presented for analysis of two-particle azimuthal angle correlation functions ob- 
tained in collisions at ultra-relativistic energies. We show that harmonic and di-jet contributions to 
these correlation functions can be reliably decomposed by two techniques to give an accurate mea- 
surement of the jet-pair distribution. Results from detailed Monte Carlo simulations are used to 
demonstrate the efficacy of these techniques in the study of possible modifications to jet topologies 
in heavy ion reactions. 
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A primary objective for studying ultra-relativistic col- 
lisions between complex nuclei is the creation and charac- 
terization of a new phase of hot and dense nuclear matter, 
where quarks and gluons are no longer confined to the 
interior of single hadrons. The existence of this Quark- 
Gluon Plasma (QGP), was conjectured based on the ob- 
servation of asymptotic freedom Q, Q and is now very 
strongly indicated by lattice QCD calculations H 0, [5|. 
The energy loss of hard scattered partons in this de- 
confined medium is predicted to be larger than that for 
cold hadronic matter H Q . Such an energy loss can re- 
sult in an apparent softening of the fragmentation func- 
tion for jets, leading to suppression in the observed yield 
of high pt hadrons and to modifications of jet topologies. 
Thus, a careful study of jet properties and yields in heavy 
ion collisions provides a sensitive probe for possible QGP 
formation H [1] and possibly for QGP properties. 

It is known that characteristic back-to-back jets of 
hadrons are copiously emitted in sufficiently energetic 
collisions. Studies of these jets provide information 
on related fra gm entation functions in a "cold" nuclear 
mcdium|9l llCl lllj . The more recent program of study at 
the Relativistic Heavy Ion Collider (RHIC) has involved 
p-p, d-Au, and Au-Au collisions. It is conjectured that 
the first produces jets in free space, the second produces 
jets that could be influenced by relatively modest initial 
state effects, and the third produces jets that could be 
strongly influenced by final-state effects such as those re- 
sulting from a reaction zone of very high energy-density 
nuclear matter. Hadron pair correlations provide a pri- 
mary tool for the study of the yield and topology of these 
jets. Therefore, much excitement revolves around the ob- 
servation, comparison, and characterization of such cor- 
relations in these three situations [l2L fl3l fTi! flEl Il6j| . 

Pair correlation studies of d-Au and Au-Au collisions 
have to deal with a relatively large "underlying event". 
For d-Au reactions, this underlying event is uncorrelated 
with the reaction plane (RP). By contrast, for Au-Au 
reactions there are rather strong correlations between 
the reaction plane and the azimuthal emission angle of 



hadrons [11111111113. This "elliptic flow" leads to 
particle-plane-particle or "harmonic" pair correlations 
that must be taken into account as one tries to iden- 
tify and characterize jet pair correlations. Flow corre- 
lations are presumed to result primarily from pressure 
gradients developed in the initial anisotropic (or almond 
shaped) overlap region of high energy-density nuclear 
matter 0, |H |2l [gl Q • The interactions within this 
overlap region may also alter jet properties, and one has 
to allow for this possibility in any analysis. 

Recent experimental studies of jet properties have fo- 
cused on the decomposition of A</> distributions (Til ITfjl 
I2H23. Here, we present two methods which exploit the 
use of correlation functions to separate jet driven pair 
correlations from harmonic correlations. The use of cor- 
relation functions circumvents the need for full azimuthal 
detector acceptance and can serve to minimize many im- 
portant systematic uncertainties (acceptance, efficiency, 
etc) which could influence the accuracy of extracted jet 
properties. A unique capability of one of the methods 
presented is the direct study of modifications to jet prop- 
erties via the extinction of harmonic correlations. The 
methods are illustrated and tested via detailed Monte 
Carlo simulations. Applications to the analysis of RHIC 
data are presented elsewhere. 

Azimuthal angle correlation functions: To em- 
phasize jet correlations, one often selects events with at 
least one high-p^ particle. This "trigger particle" is then 
paired with other "associated" particles to obtain az- 
imuthal angle pair correlations. Hereafter, we refer to the 
trigger and associated particles as A and B respectively. 
The AB correlations are termed "fixed correlations" if 
A and B are from the same px bin or "assorted correla- 
tions" if A and B are from separate pt bins. AB pairs are 
used to construct a correlation function C AB (A(f>) in az- 
imuthal angle difference Acj) = ((f) a — 4>b)] this function is 
constructed by dividing an area normalized (foreground) 
distribution of real particle pairs by an area normal- 
ized (background) distribution of mixed pairs obtained 
by pairing particles from different events. 
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The foreground distribution F is defined in terms of 
the number of coincident (or real) pairs N AB per event; 

F(Acf>) = jV AB (A0)/EiV AB (A0), (1) 

where the sum is taken over the bins in A<p. 

This gives the probability distribution for detecting an 
AB pair in a bin of A(j). The background distribution 
is defined in the same way, but by selecting A and B 
particles from separate events to form pairs; 

B(Acf>) = [N AB {A$)/VN AB {A$)] mix . (2) 

This gives the detection efficiency for an AB pair in the 
bin of Aip. The ratio of these two distributions defines 
the correlation function 

C{A4>) = F(A<f>)/B(A<t>). (3) 

This correlation function gives the efficiency corrected 
probability distribution for real coincident pairs. 

The filled circles in Fig. ^ show a typical correlation 
function generated from simulated data tuned to recent 
experimental data from RHIC [PH Wfy . There is a rel- 
atively narrow peak at 0°, a broader peak at 180°, and 
a minimum at ~ 80°. Our task is to devise an analy- 
sis method that allows one to retrieve the input jet-pair 
distribution used to generate the correlation function. 

Two source model: For Au-Au reactions at RHIC 
it has been shown that a dominant aspect of hadron pair 
correlations is collective flow of second order, i.e. ellip- 
tic flow 0, |22| ■ This gives a harmonic distribution for 
particle type A (or B) with respect to the angle of the 
reaction plane ipR. 

N A ^ A ~ i>n) oc [1 + 2v A cos2{<j )A - (4) 

and a correspondin g h armonic correlation function be- 
tween AB pairs psll^ 

C£ B {A<P) = [l + 2«2C<w2(A0)]; v 2 = (v A xv B ). (5) 

In addition to these harmonic correlations one expects 
characteristic contributions from jets. To simulate both 
source contributions we include both harmonic and di-jet 
particle emission patterns into a Monte Carlo code. The 
simulations use exponential px distributions and Poisson 
sampling for the number of jets per event, the number of 
particles per jet and the number of harmonically flow- 
ing particles per event. The average number of near 
and away-side jet particles were set to be equal, and the 
near and away-side jets were generated with an effective 
jx and kx respectively. Here, jx reflects the average 
transverse momentum projection of hadrons perpendic- 
ular to the nearside jet axis (jx = p!^ adron sin (4>j e t-had)) 
and kx reflects the acoplanarity of the near- and away- 
side jets, {k T — (PT*/VT 2 )) sm (^^et-iet))- Input pa- 
rameters were adjusted to yield correlation functions in 
rough agreement with the measured ones. An absorp- 
tion factor was introduced to study the effects of possi- 
ble suppression of the away-side jet (J a ). The axis of the 



nearside jet (J„) was also made to correlate with vary- 
ing strength to the RP via a harmonic function (as in. 
Eq. |@}) or left uncorrelated. 

It can be shown |30| that the pair correlations from 
such flow and jet sources are given by 

C AB (Ac/)) = a a [C AB ( A<t>) + J( A0)] , (6) 

where Cf^ B {A(f) is a harmonic function of effective am- 
plitude v 2 , Eq. J5j), and J(A^) is a di-jet function. No 
explicit or implicit assumption is made for the functional 
form of of J. In Fig. Q the dot-dashed and solid curves 
represent the input harmonic and jet functions respec- 
tively. 

By rearrangement of Eq. © one obtains 

J(A4>) = [C AB (A0) - a o C£ s (A0)]/a o , (7) 

Knowledge of a Q is required to evaluate J. Next we dis- 
cuss an ansatz for assigning the value of a Q . 

From measurements of the di-jet functions in d+Au 
reactions, one finds that they are very similar to the 
ones obtained in p+p reactions and can be reasonably 
described by a double Gaussian function with relatively 
narrow near- and away-side widths (<7„ and a a respec- 
tively) m l3T| . For similar px selection, Au+ Au reac- 
tions arc found to give near-side jet widths which are 
similar to those obtained in d+Au. By contrast, the 
away-side jet seems to broaden and to diminish in am- 
plitude as the centrality is increased 0|- Guided by 
these observations it is reasonable to assume, as a start- 
ing point, that the di-jet function in Au+Au reactions 
has negligible intensity at the minimum A0 m 2 n , in the 
jet function (cf. Fig. i.e. 

a Q C% B {A$ mm ) = C AB (Acj) m i n ), (8) 

which can be solved to obtain ao- Thus the condition, 
zero yield at minimum (ZYAM) , serves to fix the normal- 
ization constant a Q . 

From the two source model described above, one ob- 
tains a single particle distribution, with respect to the 
reaction plane ipR which can be fit with a harmonic. 
Its amplitude i>2, reflects an average over the particles 
from both sources. For example, if one simulates a 
flow source (90%) with V2 = 0.2 and an uncorrelated 
jet source (10%), then the observed value of will be 
0.9 x 0.2 + 0.1 x or 0.18. A finite resolution in the 
determination of ipn requires a correction factor in the 
evaluation of v-i as described in Refs. [22, EE [HI l^fjj - 
Knowledge of 1)2 determines C^ B in Eq.© The dete rmi- 
nation of the effective V2 utilizing the reaction plane as 
discussed above, requires that the reaction plane is itself 
determined by a procedure free of non-flow effects. One 
such example is a measurement which demands a large 
(pseudo)rapidity gap (^3 — 4 units) between the reac- 
tion plane and the particles correlated with it j23, H(| . 
It is expected that the V2 values so obtained are much 
less affected by jet contributions [27L l3f|. Earlier work 
has shown that V2 varies with both centrality and px- 
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Therefore, to avoid biases, one must measure v£ and 
for exactly the same data set to be analyzed for its jet 
function. To accomplish this, one also uses the reaction 
planes and their resolution corrections for exactly the 
same events of interest and determines and from 
Eq. 01 for the particle of interest. From Eq. 0one then 
obtains the effective amplitude t>2- 

With the ZYAM condition and the measured values of 
V2, one can decompose the observed correlation function 
into jet and harmonic components. This is illustrated in 
Fig. ^ The solid circles represent the simulated correla- 
tion data. The dot-dashed curve represents the effective 
harmonic component about ao- The filled triangles show 
the difference with an ao offset, which can be compared 
to the input jet function (also with an ao offset) repre- 
sented by the solid line. Note that the values of v£ and 
i>2 are obtained from fits to the azimuthal distributions 
with respect to the reaction plane. 

The method can be extended to obtain information on 
the correlation of the jet function with respect to the re- 
action plane |37l l3Sj . Two new correlation functions are 
constructed as follows: (1) the trigger particle A is con- 
strained to have its <fi angle within a certain cut-angle 



A(f> c , about the reaction plane, and the associated par- 
ticle B is unconstrained. The mixed-event background 
B(A<fi) is constructed as before in Eq. 

C m ( A<t>) = F m ( A0) /B(A4>) ■ (9) 

(2) The trigger particle is constrained to have its <f> angle 
within Acf> c of the normal to the reaction plane 

C out {A<j)) = F out (A<i>)/B{A4,). (10) 

Figure [5] shows the effect of these constraints on simu- 
lated data for pure flow with no di-jet contribution. De- 
spite a sizable dispersion of the reaction plane (At/jr = 
40°), Fig. [2| shows that the effects of these constraints 
are strong and are in opposing directions for C m (A(f>) 
versus G out (A<f)). The smooth curves in the figure show 
results derived from the analytical formulas developed in 
Ref. [37| for the relationship between (v£) in and (v£) out 
for given values of v^, the reaction plane resolution AtpR, 
and the angular selection A<j> r us ed to constrain the di- 
rection of the trigger particle [33 • 



2v£ {Acj) c ) - sin (2Acj) c ) (cos (2A4^)) + \ sin (4A0 C ) (cos (4A* H )) 
2 {A<t> c ) - 2v£ sin (2A(j) c ) (cos (2A^ R )) 

'2v£ (A(f> c ) + sin (2A0 C ) (cos (2A* 7? )) + ^ sin (4A0 C ) (cos (4A# H )) 
2 (A(j) c ) + 2v£ sin (2A<j) c ) (cos (2 A* R )) 



The harmonic amplitude Vi in Eq. [S] is replaced by 
( v 2) m i v 2) an( l by (v^Y^ivf)- I n the presence of jets, 
one can also use the ZYAM condition to decompose in- 
plane or out-of-planc correlations. This is illustrated in 
Fig. |3 where results are generated that include harmonic 
emissions of 92% of the particles and di-jet emissions of 
8%. Panel (a) shows simulated data points compared to 
a fit by Eqn. 0] for the particle distribution with respect 
to the reaction plane. This fit gives the effective value of 
V2 to be used in the analysis of the two-particle correla- 
tion functions. Panels (b), (c), and (d) show simulated 
data for an inclusive, in-plane and out-of-plane correla- 
tion function (filled circles) respectively. The harmonic 
correlations in each figure (dashed-dot curve) are calcu- 
lated analytically from Eqn. 0] using i>2 from the fit in 
panel (a) along with (t>^) OMi and (v^) 1 " from Eqs. ITTI 
and 1 121 The ZYAM condition is used to fix the constant 
ao (in each case) as described above. Then the extracted 
points for jet-pair distributions are determined by differ- 
ence via Eq.[7| Comparison to the input jet-pair distribu- 
tion (solid curves) shows an excellent agreement, thereby 
giving confidence in the decomposition procedure. It is 
noteworthy that a similar analysis performed for the case 



in which di-jets were made to correlate with the reaction 
plane, gave equally good results. 

As indicated earlier, this method of decomposition 
makes no assumption about the functional form of the 
jet function. From the shape of the jet functions ob- 
served in p+p and d+Au reactions, one is led to expect 
a double Gaussian jet function peaking at and 180 deg. 
However, if the reaction creates a strongly interacting 
medium, then that medium could distort the jets, espe- 
cially on the away-side. Such distortions have been re- 
cently attributed to to the influence of collective flow [3^ 
and the generation of shock waves around partons prop- 
agating through the medium j4j}. We have performed 
several simulations in which strongly distorted away-side 
jets were introduced. In each case the decomposition 
method retrieves the input jet function in detail, con- 
firming that the decomposition procedure is robust even 
for unusual di-jet distributions. 

Extinction of Harmonic Correlations: Figure |3| 
shows pair correlations with strong harmonic components 
C^ s (A0) and relatively weak jet components J(A</>). 
Only after subtraction of C^ B (A</>) is the true shape re- 
vealed for J(A0). It would be desirable to select a data 
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set that contained only those correlations due to the jets; 
Fig. |21 and Eqs. [S] and ^2 point the way. As discussed 
above for Fig. the out-of-plane constraint on particle 
A changes the phase of the harmonic correlation to favor 
pairs at A<f> = 90°. By using Eq. ^]one can select the 
particular cut angle A(j) c = 4> xt such that (v^) OM * = 0. 
For this case the value of V2 is driven to zero in Eq. [5] 
and the harmonic correlations arc extinguished. 

In Fig. 0] we demonstrate the extinguishing of har- 
monic correlations. The filled circles show the inclusive 
A<j) distribution obtained for a simulation with an input 
i>2 ~ 0.16. The filled squares show the out-of-plane Acf> 
distribution after the cut angle <f> c , is set to the extinction 
value, ie. </> c = (f> x t. The latter distribution is flat and 
demonstrates that this technique gives a good method to 
determine the jet correlations directly from a data set, 
without the blurring effect from harmonic correlations 
mediated by the reaction plane. The latter technique 
provides clear advantages for future jet studies at RHIC. 

Jet pair fraction and conditional yield: The sum 
of the jet-pair function J over A(f> can be related to the 
average fraction of jet-correlated particle pairs per event, 
hereafter termed jet pair fraction JPF, as follows, 

JPF = Ea J(Acf>)/T,C(A4>). (13) 

The hatched area in Fig. ^ illustrates this sum. This 
value of JPF represents the fraction of all mid-rapidity 
AB pairs that are correlated by jets i.e. those correlated 
by effects over and above those mediated by the reaction 
plane. 

One may ask what is the event-averaged number of jet- 
associated particles per trigger particle. This quantity 
is commonly referred to in jet-analyses as per-trigger or 
conditional yield. The conditional yield can be calculated 
from the jet pair fraction as follows. 

For a given data set, one records the detected number 
of AB pairs per event n AB , the detected number of A 
particles per event n A , and the detected number of B 
particles per event n B . The ratio R p of the pair rates to 
the product of the single particle rates 

Rp = n AB /(n A x n B ) (14) 

carries information about pair-production over and above 
the combinatorial background. These detected rates 
should be corrected for efficiencies to get the true rates 
for particle-pairs and singles at all <f> but with the relevant 
77 acceptance. Assuming that the pair-efficiency approx- 
imately factorizes into the single particle efficiencies, we 
can denote the true values per event by; 

n A = (n A )(Ef A ) (15) 
n B = (n B )(Ef B ) (16) 
n AB « (n AB )(Ef A )(Ef B ) (17) 

i.e. n AB /nfnf = n AB /n A n B . 



Now one writes for the true efficiency corrected condi- 
tional yield TCY 

TCY = JPF(n AB /{n A x nf )) x nf (18) 
= JPF x R p x nf (19) 

From the procedure outlined above one gets JPF, and 
one has R p from the detected values n AB , n A and n B 
for each data set. From the literature (see for example 
Refs. ^1E3) one has previously published pr spectra for 
many particles corrected for azimuthal efficiency. These 
can be integrated over the px intervals for particle B to 
get nf . 

The evaluation of systematic errors in this procedure 
will center on the determination of the jet pair fraction 
(shaded area in Fig. ^ for example). This is because the 
jet-pair distribution could have a finite (albeit small) in- 
tensity at A(f> m i n , in certain situations. If this is indeed 
the case, the extracted jet shape is essentially unaffected, 
but the yields extracted under the ZYAM assumption 
would constitute lower limit yields. On the other hand, 
one can adopt an appropriate functional form for the 
shape of the jet function to make an estimate of the pos- 
sible jet yield at A</> m j„, and hence, assign a systematic 
error via iteration. 

In the discussion of jet yields we used Eq. H17f) 

n AB = n AB (Ef A )(Ef B ). 

In fact, we note that AB pairs that are correlated in A<j> 
may well have a somewhat modified efficiency, designated 
below by M,. 

n AB = nf B Ef A Ef B (M l ) (20) 

After determination of the harmonic amplitude one may 
get the Mh factor by Monte Carlo simulation techniques. 
Similarly one can get modification factors M. n j and Wl a j 
for the near-side jets and the away-side jets respectively. 

Correlations between the multiplicities n A and n B 
make the mean number of combinatoric pairs n AB nb 
greater than the product n A n B . One can write 

ni B Jn A n B t = l, £ jL 1. (21) 

However the quantity £ cancels out in the product 
(JPF)(n^ B ) in the present analysis. 

Summary: A Monte Carlo reaction code has been used 
to generate simulated data sets for hadron pair correla- 
tions tuned to the general behavior from Au + Au reac- 
tions at RHIC. These data sets contain pair correlations 
from both elliptic flow and from jets. A method is pre- 
sented for the analysis of such data to retrieve the jet 
pair fractions, conditional yields, and widths for both 
near and away side jets. A novel technique is also pre- 
sented for direct quenching of harmonic correlations to 
reveal jet correlations. 
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FIG. 1: Correlation function from simulated data (filled cir- 
cles). The dotted curve shows the input harmonic component. 
The solid curve and full squares show input and output jet 
pair distributions referenced to an offset. 




FIG. 2: Simulated data, closed and open circles, for inplane 
C m and out of plane C out correlation functions as described 
in the text with A<f> c = 45° and Aipa = 40°. Only harmonic 
correlations are present as shown by the dotted curve with 
amplitude v% = ■ 
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FIG. 3: (a) Simulated data points and smooth curve for 
the correlation with respect to the reaction plane, (b) Simu- 
lated data (filled circles) for an inclusive correlation function 
including harmonic and jet correlations . The dotted curve 
indicates the harmonic term. The solid curve and full squares 
show input and output jet distributions referenced to a . (c) 
Same as (b) but for an in-plane correlation function Eq.(9). 
(d) Same as (b) but for an out-of-plane correlation function 
Eq.(10). 
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FIG. 4: Simulated correlation function for unconstrained 
particles (filled circles) and for a trigger particle constrained 
within the cut angle cj> c = cj> x t perpendicular to the reaction 
plane (filled squares), see text. The results are for a pure 
harmonic simulation with V2 ~ 0.16. 



